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Abstract. 

We model the OH megamaser emission from the luminous infrared galaxy IIIZw35 as arising from a narrow rotating starburst 
ring of radius 22 pc enclosing a mass of 7x10* Mg. We show how both the compact and apparently diffuse maser emission 
from this ring can arise from a single phase of unsaturated maser clouds amplifying background radio continuum. The masing 
clouds are estimated to have a diameter of < 0.7pc and internal velocity dispersion of ~20kms"'. We find that the clouds are 
neither self-gravitating nor pressure confined but are freely expanding. Their dispersal lifetimes may set the vertical thickness 
of the ring. For an estimated internal density of 3 x 10^ cm"^, cloud masses are of order 24Mo. The observed spectral features 
and velocity gradients indicate that the clouds must be outflowing and escaping the nucleus. The cloud mass outflow rate 
is estimated to be O.SMoyr"', while the star formation rate is ~19MQyr"'. Associated ionised gas, possibly generated from 
dissipated clouds, provides free-free absorption along the source axis, explaining the observed East- West asymmetries. We 
show that the dumpiness of a maser medium can have a dramatic effect on what is observed even in a relatively low gain OH 
megamaser. Specifically, in IIlZw35 our clumpy maser model naturally explains the large line to continuum ratios, the large 
1667MHz: 1665MHz line ratios and the wide velocity dispersions seen in the compact maser spots. Other astrophysical masers 
showing both compact and apparently diffuse emission might be explained by similar clumpy structures. 

Key words. Galaxies: starburst - Galaxies: individual: lllZw35 - Masers 



1 . Introduction 

Extra-galactic OH megamaser emission is generally associated 
with compact (< 100 pc scale) starburst or AGN activity in 
the centres of IR luminous galaxies. Observations of this maser 
emission provides a unique method of studying the structure 
and kinematics of galactic nuclei at parsec resolution without 
dust obscuration effects. Measurements of velocity gradients 
and line widths (e.g. Pihlstrom et al. (2001), hereafter POl) al- 
ready provide important constraints on stellar mass densities 
and turbulent velocities in IR luminous galaxies. Such obser- 
vations may also trace large scale obscuring tori in composite 
AGN/starburst sources (Klockner et al. 2003). Potentially OH 
megamasers can also tell us about the size, density and tem- 
perature of molecular clouds in the central ISM of starburst 
galaxies and AGN. However, to accomplish these goals a better 
understanding of the OH megamaser phenomenon is required. 

Recent MERLIN and VLBl observations of OH mega- 
masers have suggested that the standard model, developed in 
the 1980s (Baan 1989), might need modification (see Lonsdale 
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(2002) for a review). In this standard model OH maser emis- 
sion is generated by low gain (|r| ^ 2) unsaturated amplifica- 
tion of background continuum by a foreground OH amplifying 
medium. Although this medium is comprised of discrete OH 
clouds it is implicitly assumed that there are many such clouds 
and they individually have very low gains. These clouds there- 
fore form an effective gas in which statistical fluctuations in 
cloud number between different lines of sight are unimportant 
and so the maser opacity varies slowly across the source. Given 
these assumptions the amplifying medium is often described as 
a diffuse screen. 

Consistent with the standard model, early VLA and 
MERLIN observations showed that OH maser and continuum 
emission overlapped. Since both continuum and maser emis- 
sion were assumed to be smooth on scales ^ 100 pc VLBI 
observations were not expected to reveal anything interest- 
ing. However, when VLBI observations were finally made 
of Arp220, both compact continuum (Diamond et al. (1989), 
Smith et al. (1998)) and compact OH maser emission were 
detected (Lonsdale et al. (1994), Lonsdale et al. (1998) here- 
after L98). Remarkably, the bright maser spots in Arp220 
were not coincident with the continuum spots, and some maser 
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spots displayed extreme line-to-continuum ratios (i.e. > 800 
in L98). These observations were clearly inconsistent with dif- 
fuse screen models. Slightly less extreme compact maser emis- 
sion was subsequently detected in other sources (Trotter et al. 
(1997), hereafter T97, Diamond et al. (1999), hereafter D99, 
Klockner et al. (2003), Klockner & Baan (2004)). These same 
sources also contain diffuse maser components which account 
for between 50% and 90% of the total maser flux density. It 
has been suggested that the compact masers occur in saturated, 
perhaps collisionally pumped regions, while the diffuse compo- 
nent comes from an extended unsaturated, radiatively pumped 
screen fully consistent with the standard model (L98, D99). 

One of the clearest cases of an OH megamaser showing 
both compact and diffuse maser emission is IIIZw35. Two 
groups of compact masers were detected in VLBI observations 
(T98, D99), but only about half of the total OH maser emis- 
sion seen on MERLIN scales (Montgomery & Cohen 1992) 
was recovered in the VLBI maps. To determine the location of 
this missing component POl conducted EVNh-MERLIN obser- 
vations. These observations revealed perhaps the best example 
of a rotating OH maser ring yet found. The previously known 
compact masers lie at the tangents of this ring. Such a location 
is hard to understand if compact and diffuse masers are gen- 
erated by different physical mechanisms. Instead, a geometri- 
cal origin for the compact masers is suggested. POl proposed a 
mechanism based on a single phase of OH masing small clouds 
(~1 pc). At the ring tangents multiple overlaps between clouds 
in space and velocity are likely due to the increased path length 
through the ring. These multiple cloud overlaps give rise to the 
bright compact maser features. Elsewhere, at the front and back 
of the ring where there are few such cloud overlaps, the emis- 
sion consists of many weak maser spots. These spots are too 
weak to be detected individually in high resolution observa- 
tions but in low resolution observations they are averaged to- 
gether and give rise to an apparently diffuse emission. 

In this paper we investigate in greater depth the clumpy 
ring model proposed in POl. Here we fully consider the spec- 
tral properties of the clouds and their velocities are treated in 
a more realistic way. We also model the continuum emission 
in a geometrically and physically consistent manner. Using nu- 
merical simulations we demonstrate that most of the available 
OH maser observations of IIIZw35 can be explained using an 
improved version of the POl model. Most of the input parame- 
ters for these simulations are constrained directly by the obser- 
vations. Our modelling also illuminates general properties of 
maser ring geometries and clumpy maser media. 

The detailed plan of this paper is as follows. Section 2 
summarizes the relevant observations of IIIZw35. In Section 
3 we present the clumpy ring model and derive its parame- 
ters. Section 4 describes the numerical simulations and com- 
pares the results of the simulations with observations. Section 5 
presents a general discussion on how the observed maser char- 
acteristics can be explained by a clumpy medium. In Section 
6 the derived physical properties of the clouds and the galactic 
nucleus are discussed. Finally in Section 7 we draw conclu- 
sions and suggest future work. 



2. Source properties and observations 

IIIZw35 is a Luminous Infra-red Galaxy (Lfir = IO'^^Lq) 
at a distance of 110 Mpc (assuming Hq = 75 kms^'Mpc"', 
(Chapman et al. (1990), hereafter C90)) which gives a linear 
scale of 0.5 pc mas"'. From optical spectroscopy it has been 
classified as a borderline LINER/Seyfert 2 (C90). The galaxy 
lies on the well known radio-FIR correlation and both the ra- 
dio continuum and IR emission are consistent with being pow- 
ered by a starburst with a star formation rate (SFR) of ~I9 
Moyr-i (POl). Both the 1667 MHz and weaker 1665 MHz 
megamaser emission were discovered by Staveley-Smith et al. 
(1987). The ratio of peak opacities in the 1667:1665 MHz spec- 
tra is ~9 (Staveley-Smith et al. (1987) and Mirabel & Sanders 
(1987)), which is at the higher end of observations amongst OH 
megamaser galaxies (Henkel & Wilson 1990). From single dish 
CO observations, Mirabel & Sanders (1987) estimate an over- 
all molecular gas mass of 6.5x10^ Mq. The dynamical mass 
within a radius of 22 pc is estimated by POl to be ~7xl0'' Mq, 
implying that most of the molecular gas is located in a more 
extended structure. 

Intermediate resolution maps of the OH maser emission 
have been produced by POl by combining EVN and MERLIN 
datasets and are reproduced in Figure 1. The greyscale in 
Figure la shows that, at this resolution, the maser emission 
has a clear void near its centre and two bright regions to the 
North and South, with an estimated Line-to-Continuum Ratios 
(LCR) of 47 ± 13 and 73 + 8 respectively. Two bridges of dif- 
fuse emission connecting the bright regions can also be seen. 
Although the one to the West is brighter, it has the smaller LCR; 
the western bridge has LCR = 9 ± 2, the eastern bridge has an 
LCR = 14 + 4. Figure lb shows the North-South velocity gra- 
dient obtained from the maser moment map. This gradient and 
the maser brightness distribution, are consistent with a rotat- 
ing ring whose axis is inclined by ; ~ 60° from the direction 
of the observer The bright North and South regions are in- 
terpreted as arising from the tangents of the ring, where path 
lengths are longer, while the bridge emission originates at the 
front and back sides of the ring (see Figure 2). At VLBI res- 
olution (T97, D99), the diffuse bridges are not detected while 
the two bright regions break up into complexes of bright maser 
spots. In one of the spots the inferred LCR is > 500 (D99). 
In general the spots are unresolved, thus having sizes of less 
than 0.7 pc (D99). Spot spectral widths are measured to be 30- 
50 km s '(D99). Both, T97 and D99 find that the spot velocity 
centroids trace a roughly East- West gradient within the north- 
ern region. Approximately half of the velocity integrated line 
flux density seen in Figure I a within the northern and southern 
regions comes from the compact spots, the rest is not detected 
at either VLBA (T97) or global VLBI (D99) resolution. 

In addition to the OH maser emission POl also detected 
continuum emission in their EVN+MERLIN data (see contours 
in Figure la). In contrast at global VLBI resolution only upper 
limits to the continuum have been found (T97; D99), implying 
that most of the continuum emission is smoothly distributed. 
The continuum emission seen at EVN+MERLIN resolution, 
just like the OH emission, is stronger on the western side of 
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Fig. 1. Images reproduced from POl showing observations of OH maser and continuum emission at EVN-l-MERLIN resolution a) Greyscale 
sliows the velocity integrated OH emission at a resolution of 34x29 mas. The contours show continuum emission at the same resolution. 
Contour levels are at -1,1,2, and 4 times the 3cr noise of 0.42mJy/beam. b) Corresponding OH maser velocity centroid field. The greyscale is 
between 8224 (lightest grey) and 8300 km s"' (black). Contours are from 8230 km s"' increasing at 5km s"' intervals up to 8290 kms"'. 



the source and is therefore asymmetrically distributed about the 
major axis of the inferred OH maser ring. 

3. Clumpy ring model 

3.1. Overall Geometry 

In this paper we fit the IIIZw35 observations by an inclined 
axisymmetric model in which both OH clouds and continuum 
emission coexist within circumnuclear rings (see Figure 2). To 
reconcile such a symmetric geometry with the observed East- 
West asymmetry in both line and continuum emission (see 
Figure la) our model also includes a bi-cone of free-free ab- 
sorption which covers the eastern side of the source. This ob- 
scuration defines the maser ring orientation requiring the east- 
em side to be the most distant. Note that the existence of a 
free-free absorbing component is supported by independent 
evidence. C90 found that the integrated radio spectrum has 
a turnover at ~1 GHz. Given the size of the radio emission 
this cannot realistically be due to synchrotron self-absorption, 
because enormous departures from particle/field equipartition 
would be needed. Physically the free-free absorbing cone could 
be the base of an outflowing superwind such is often observed 
in energetic starbursts (Heckman 2003). 

As noted in §2 although the total maser emission is weaker 
on the eastern side the LCR is about factor of 2 larger there 
than on the western side. This difference is explained if more 
of the seed continuum emission is background to the OH ring 
on the east side. This naturally occurs if most of the bright con- 
tinuum emission comes from a larger radius than the OH mas- 
ing gas (see Figure 2a). This is corroborated by the fact that 
the observed continuum radio source is larger than the mega- 
maser source (see Figure la). As shown in §4.3 this geometry 



also predicts that the brightest masers do not occur exactly at 
the maser ring tangents but slightly to the east of the tangent 
points, just as found by the observations of T97 and D99. 

The detailed properties of the different components of the 
model are estimated in the following subsections and sum- 
marised in Table 1. These parameters are used in the Monte- 
Carlo simulations described in §4. 

3.2. Free-Free absorption Bicones 

From the observed ratio of the continuum emission brightness 
on the East and West sides of the source, the line of sight opac- 
ity through the cone to the eastern bridge is estimated to be 
~0.7 at the observing frequency. Since all of the eastern bridge 
emission is weak compared to the western bridge the opening 
angle of the bicone must be significant. We therefore assume 
a bicone in which the opacity per unit length decreases with 
angle from the symmetry axis as a gaussian with FWHM 60°. 

3.3. Ring dimensions 

The dimensions of the region containing the OH masing clouds 
are shown in detail in the right panel of Figure 2. This thin ring 
structure has inner radius Ri„ and an outer radius /?out = ^in -i- 
AR. POl find that Rin ~ 22 pc. From the North-South extent 
of the region covered by bright maser spots in the VLBA and 
global VLBI maps (T97 and D99), A/? is ~ 3 pc. From its East- 
West extent near the ring tangents we infer a ring height of AH 
~ 6pc. The aspect ratio of the projected ring shows that the ring 
axis is inclined at angle / ~ 60° from the observer's direction 
(POl). We find that an axisymmetric continuum emissivity pro- 
file which peaks at radius Rc = 37 pc best fits the continuum 
observations (see Figure la and §4.2). We also find that the con- 
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Fig. 2. a) Sketch of the proposed source geometry. The inner darlc grey ring represents the region where OH masing clouds are confined. The 
outer light grey ring depicts an isosurface of the smoothly distributed continuum emissivity. Note that although the radius of peak continuum 
emissivity lies outside the OH maser ring, some continuum emission interpenetrates and even lies within the OH ring (not shown in this 
figure). This geometry explains the large line to continuum ratio on the eastern side of the source because here the majority of the continuum 
is background to the OH. In contrast on the western side only the smaller fraction of the continuum which interpenetrates and lies inside the 
OH maser zone is available as a source of seed photons. In order to explain the relative weakness of the absolute brightness of both line and 
continuum on the eastern side of the source (see Figure 1) the model also includes a region of free-free absorption within a bicone which covers 
the eastern side of the source, b): Detailed representation of the OH maser ring indicating positions and dimensions referred to in the main text. 
Arrows indicate components of cloud rotation around the ring (vrot) and outflowing from the ring midplane (I'z). See §3.6 for more details. 



tinuum emissivity function must be relatively wide in radius so 
that not all of the continuum lies outside the OH maser region. 

3.4. Need for a single-phase clumpy medium 

Having constrained the geometry of the maser emitting zone 
we now argue that it must consist of distinct clouds. Maser am- 
plification of background continuum gives LCR = e'^, where t 
is the maser optical depth. When the amplification is due to 
a smooth distribution of uniform unsaturated maser gas, r is 
proportional to the path-length through the region. The ring 
geometry derived in the previous sections implies a factor of 
8 difference between the maximum path-lengths in regions N 
and E (see Figure 2) such a ratio in path length implies LCRn 
= (LCRe)^. In contrast, the ratios observed at EVNh-MERLIN 
resolution are LCRe ~ 12 and LCRn < 75, a conflict of at 
least four orders of magnitude with this relation. The problem 
becomes even worse if saturation affects the bright features be- 
cause the optical depth of a saturated maser is always larger 
than that of an unsaturated maser of the same length. 

In contrast with a smooth maser screen, a collection of in- 
dependent clouds involves a statistical effect, predicting a vari- 
ation in LCR with position around the ring which can match the 
EVNh-MERLIN observations. As described in the appendix, a 
remarkable property of clumpy media is that the effective opac- 
ity is not linearly proportional to the path length, instead, it falls 
below the linear prediction. This result gives a compelling ar- 
gument for the OH masers occurring in clouds. This argument 



is separate and additional to the one proposed in POl which re- 
quired clouds to explain the brightest spots by multiple cloud 
overlaps. 

3.5. Cloud Properties 

Based on the high resolution observations of the brightest 
maser spots in D99 maser spot diameters are estimated to be 
< 0.7 pc. In general the spot size associated to a high gain 
maser is smaller than the size of its originating maser cloud 
(by a factor of ~ -\/t). However since the maser opacity of the 
clouds employed in our model is about one (see below) the size 
of the maser produced by a single cloud is about the same size 
of the cloud. In spots caused by multiple cloud overlaps the ap- 
parent reduction in size caused by increased opacity is counter- 
acted by the widening due to spatial misalignments. The cloud 
size is therefore estimated to be the same as that set by spot 
size, < 0.7pc. A lower limit for the cloud size is set from physi- 
cal arguments (see §6.1). For our Monte-Carlo simulations (§4) 
we adopt the observed upper limit as the actual cloud size but 
argue that our results are only weakly dependent on the exact 
size chosen. 

From the observed velocity widths of the bright maser spots 
(D99) cloud internal velocity dispersions of order 20 kms ' 
are implied. Again a similar argument applies in velocity as in 
space to explain why we expect observed and true cloud veloc- 
ity widths to be similar. Values for the cloud number density 
and opacity were found from the Monte-Carlo simulations. As 
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Table 1. Model Parameters 



Maser and continuum ring properties: 

Maser ring inner radius i?in'' 22 pc 

Maser ring radial thickness Ai?" 3 pc 

Maser ring height AW 6 pc 

Maser ring Rotation velocity v.ot'' 57 km s"' 

Continuum peak emissivity radius Rj' 37 pc 

Inclination angle 60° 



Free-free absorption bicone properties: 
Opacity T(f 0.7 
Opening angle (FWHM) 60° 



Cloud properties: 

Cloud size a*-' <0.7 pc 

Internal velocity width 6\^-'^ 20 km s"' 

Random velocity dispersion Av'^ '' 60 km s"' 

Outflow velocity V; 60kms"' 

Outflow velocity dispersion Av-^ 30 km s"' 

Number of clouds along radius A'o 1.8 

Cloud volume filling factor / < 0.08 

Cloud maser optical depth T(. 1.5 



Refs: ''T97, ''POI, =^^099 

The dispersions represent the FWHM of a 
gaussian probability distribution. 



described in §4.1, the initial starting point for this search was 
guided by the observed LCR distribution around the ring (see 
§3.4) and the theory developed in the appendix. The best fit was 
obtained assuming a cloud opacity of Tc = 1 .5 and a volume 
number density such that the number of clouds along a radius 
in the disk plane at any velocity was No - 1-8. For a cloud size 
of fl < 0.7pc this implies, averaged over the OH maser ring, a 
cloud volume filling factor / < 0.08. As we argue in §4.1 our 
fits to the data are expected to remain similar as a function of 
a as long as the cloud number density is adjusted to maintain 
a constant A^,,. In this case the cloud volume filling factor will 
scale with cloud size as / = 0.08(a/0.7pc). 

3.6. Velocity field 

Observations show (see Figure lb) a large scale North-South 
gradient in maser velocity centroid consistent with ring ro- 
tation. Note however that the velocity dispersions within the 
bright northern and southern tangent regions are comparable to 
the rotation velocity so that the integrated spectra from the two 
regions overlap (see POI and T97). At higher resolution global 
VLB! observations show within the northern tangent region a 
clear velocity gradient which is almost East-West. This small 
scale gradient is almost parallel to the projected ring axis and 
almost at right angles to the gradient found on large scales. 

The observed large scale gradient can readily be fitted if 
the clouds are assumed to have an orbital velocity component 
(see Figure 2b) of v,ot = 57 kms"' (POI). Somewhat harder to 
explain is the gradient seen within the northern tangent region. 
Part of this gradient can be explained by the ring rotation. The 



brightest masers are not found exactly at the OH maser ring tan- 
gents but slightly to the East where there is more background 
continuum (see Figure 2a) and the orbital velocity gives a small 
velocity gradient along the ring. This mechanism cannot how- 
ever explain the magnitude of the observed gradient. Instead 
we propose that clouds, in addition to their rotation velocity, 
have a comparable outflow velocity (v^, see Figure 2b) which 
is parallel to the ring axis and directed away from the disk mid- 
plane. To see how this gives the observed gradient consider a 
cloud near a tangent region in Figure 2b. If such a cloud Ues 
above the ring midplane (near side), the cloud will be located 
to the East of the tangent point and its projected outflow veloc- 
ity will be blueshifted. On the other hand if it lies below the 
ring midplane (far side) the cloud will be located to the West 
of the tangent point and its projected outflow velocity will be 
redshifted. This mechanism can produce an apparent East- West 
velocity difference of nearly twice the projected outflow veloc- 
ity within a small distance comparable to the maser ring height. 

It is interesting that at the northern tangent the cloud out- 
flow mechanism gives a velocity gradient in the same direction 
as that caused by the rotation mechanism, thus reinforcing the 
gradient. In contrast in the South the gradients from the two 
mechanisms have opposite directions and will partly cancel. 
This may explain why an East- West gradient is only seen in the 
North. Alternatively the difference may mainly be due to the 
statistical nature of our clumpy maser model. We have found 
using Monte-Carlo experiments (see §4) and assuming an out- 
flow velocity of Vz - 60 km s"' and dispersion Avz - 30 km s ' 
that small scale gradients of the required amplitude often arise 
on the tangent regions, with the northern region being favoured. 

While the outflow velocity component contributes signifi- 
cantly to the high velocity dispersion seen at EVNh-MERLIN 
resolution we also find from our Monte-Carlo modeling that 
an additional cloud random velocity component of dispersion 
Av = 60 kms ' is required. In particular this random veloc- 
ity component is needed to fit the observed very wide veloc- 
ity dispersion of the northern and southern tangent regions at 
EVNh-MERLIN resolution. 

4. Modelling and results 

In this section we calculate in detail the emission from the 
clumpy maser ring using Monte Carlo simulations. In §4.1 we 
present the details of the model. In §4.2 we explain how we 
compute our synthetic OH maser cubes and continuum images. 
Finally in §4.3 we compare our results with observations. 

4.1. Model Details 

Our numerical model contained three components (see §3.1): a 
high brightness temperature radio continuum emission to pro- 
vide seed photons for maser amplification, a free-free absorb- 
ing cone and a randomly distributed population of OH masing 
clouds. The continuum component was modelled as a radially 
smooth axisymmetric emissivity function. This function was 
chosen to have peak emissivity at radius /?c (see Table 1) and 
its shape was chosen so that it approximately matched the POI 
continuum image (see Figure la). The bi-cone of free-free ab- 
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sorption (see Figure 2a) was assumed to have a density obey- 
ing a Gaussian distribution of FWHM=60° around the ring axis 
producing a total opacity of t = 0.7 toward the eastern bridge 
region (see §3.2). 

The maser clouds were assumed to be identical and spher- 
ically symmetric with a number density of maser molecules 
described by a three dimensional Gaussian distribution with 
a FWHM of 0.7pc (see §3.5). For computational reasons this 
density profile was cutoff beyond a radius of Ipc. The corre- 
sponding cloud internal velocity profile was assumed to be a 
Gaussian with FWHM 6v = 20 km s^' . 

As discussed in §6.1 the clouds could be physically smaller 
than the above FWHM of 0.7pc and still give the required 
maser opacity. Fortunately we find that when the clouds are 
small enough to be unresolved by the highest resolution inter- 
ferometer beam the model results are only very weakly depen- 
dant on the actual cloud size used. As described in the appendix 
the critical parameters for determining structure are the cloud 
opacity and mean number of clouds per line of sight .^V. If 
a smaller cloud size was used and consequently the volume 
density of clouds increased to maintain the same A^, then the 
results to first order would be the same. The only explicit de- 
pendence on cloud size is via the inverse Poisson function (see 
appendix) which is an extremely weak function of the num- 
ber of independent lines of sight (M); a quantity which in turn 
is inversely proportional to the cloud area. The effect of using 
clouds much smaller than used in our simulations would be to 
make the peak line to continuum and 1667; 1665MHz line ratio 
somewhat larger, but as we shall see (see §4.3) these ratio limits 
are well fitted even when using the largest possible cloud size. 

As described in §3.6 the cloud bulk motions comprise or- 
bital, outflow and random components. To implement these 
motions in our simulations we first found for each cloud six or- 
bital parameters for a random orbit around the gravitational po- 
tential induced by the 7xlO^Mo enclosed mass. The maximum 
allowed orbit inclination was set to the ratio AH/2Ri„. Clouds 
were also constrained to lie within an annulus with inner and 
outer radii Ri„ and ^out=^in+ AT? (see Table 1). The outflow 
from the ring plane (see §3.6) was simulated by adding a veloc- 
ity component parallel to the ring axis to each cloud. After ex- 
perimenting it was found that assuming a mean = 60km s ' 
with dispersion Av- = 30km s"' directed away from the ring 
plane, plus a random 3D velocity of 60km s 'gave acceptable 
results. This velocity field produced realisations which showed 
velocity gradients within the northern tangent region and could 
also match the observed velocity dispersion. 

As described in §3.5, the volume number density of clouds 
Nc and their opacity Tc were left as adjustable parameters 
within the Monte-Carlo simulation. To guide the search for 
these parameters we utiUsed the formalism developed in the ap- 
pendix. A system of two equations was set up using the expres- 
sions for the beam-averaged gain over the northern and eastern 
regions. The system was solved in the least-squares sense sub- 
ject to a lower-bound constraint established by the observed 
peak gain at global VLBl resolution. The optimum solution 
found had - 1.5 and a cloud number density in the mid- 
plane of Nc = 0.2 pc^^*. For this density the number of clouds 
intersected at any velocity by a radial equatorial ray is A^o=l-8. 



Table 2. Observed and modeled quantities 





Observed 


Model 


Single dish: 






Velocity range (km s ) 


2/0 ±20 


~280 


EVN+MERLIN: 






T"North 


3.85+0.34 


3.58 


^South 


430+0.12 


3.80 


TEast 


2.66±0.33 


1.20 




2.26±0.21 


0.69 


Global VLBI: 








>6.21 


7.28 


FWHM of peak spot (km s"' ) 


-30 


30 



"Velocity width of the 1667 feature at 10% of the peak 
flux density. 

•"The effective opacities are TNom, = In(LCRNorih), etc. 



The required number density was achieved in our simulations 
by using 1200 clouds. 

4.2. Model Computation 

Maser emission was calculated by ray tracing for each line 
of sight (LOS) and velocity taking into account the amplifi- 
cation of continuum seed photons by the interspersed maser 
clouds. Synthetic spectral line data cubes and continuum im- 
ages were produced and convolved with appropriate beams for 
comparison with observations. The simulations were run in a 
dual Pentium 111 1 .5GHz computer with 2GB of RAM running 
Linux. A realisation using a synthetic cube of size 20 1^* with 
1200 clouds took approximately 20 min. An important aspect 
of our implementation is that for each realisation, the cloud in- 
formation is stored as a list of real numbers (and not gridded), 
so it is relatively easy and efficient to recompute the expected 
emission at different resolutions by specifying different grid 
sizes. 

4.3. Results and comparison with observations 

A set of 100 simulations were run using the parameters given 
in Table 1. The maser and continuum emission from a typical 
realisation are shown in Figure 3 in a way that can be directly 
compared to the observations in Figure 1 . Relevant numerical 
results for this realisation compared to observations are given in 
Table 2. The model continuum emission shown by the contours 
in Figure 3a matches well the high brightness regions of the ob- 
served continuum (compare to contours in Figure la). It is clear 
that our continuum model does not include the very extended 
emission present in the observations. However we consider that 
this model is sufficiently accurate for our main purpose of cal- 
culating synthetic OH maser emission cubes. 

The greyscale in Figure 3a shows the integrated OH maser 
emission from our model and it can be seen to be very similar 
to the observations (compare to Figure la). Figure 3b shows 
the model velocity field obtained from the moment map of the 
OH emission. Again, its overall structure agrees quite well with 
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Fig. 3. a) Model continuum emission superimposed on the velocity integrated OH maser emission. The contours, greyscale and the resolution 
are equivalent to those in Figure la. b) Corresponding modelled OH maser velocity centroid field. Greyscale is between -60 and 60 km s"' around 
the systemic velocity. Contours are from -30 km s"' and increasing by 5 km s"' up to 30km s"' . Compare with Figure lb. 



observations (compare to Figure lb). There is a North-South 
asymmetry in the model velocity field which in this realisation 
is partly caused by the presence of a single bright maser fea- 
ture in the South. The outflow velocity given to the clouds also 
contributes to break the North-South symmetry (see §3.6). The 
observations (Figure lb) show a similar but larger asymmetry 
between North and South. This might be explained by the fact 
that in the observations the south contains an even more domi- 
nant single spot than is the case in the model. Additionally note 
that the model moment map is produced using all the maser 
emission from the high resolution cube whereas in the obser- 
vations there may be a resolved-out component not included in 
the average. 

In order to show the range of structures produced in the 
Monte-Carlo simulations Figure 4 displays three selected reali- 
sations of the model and the average over the entire set. The top 
two rows show velocity-integrated images at MERLINh-EVN 
and global VLBI resolutions, the bottom row shows the corre- 
sponding integrated spectra. As expected, the simulations show 
bright emission at the tangent regions and diffuse emission in 
between. The intensity is lower in the eastern bridge due to the 
effect of the free-free absorbing cones, yet the LCR is higher 
there, reflecting the relative locations of the maser clouds and 
radio continuum. From Table 2 we see that the values obtained 
for the EVNh-MERLIN optical depths in the bridges are some- 
what less than observed. This may be due to the fact that the 
model uses only one type of spherical cloud representing a 
whole population that certainly has various sizes and optical 
depths. The bridge regions have small effective optical depths 
and can be expected to be more strongly affected by fluctu- 
ations, which are amplified exponentially, around the mean 
cloud optical depth. At global VLBI resolution (Fig 4 mid- 
dle row) the bright tangent emission breaks up into numerous 
compact maser spots due to multiple cloud overlaps. With the 



greyscale level chosen there are more such spots on the western 
side of the projected ring major axis, this is caused by the effect 
of the free-free absorbing bicone. 

The single dish model spectra (Fig 4 bottom row) agree 
in shape and velocity width with the Staveley-Smith et al. 
(1987) observations. The spectra averaged over the northern 
and southern tangent regions (dashed and dotted lines) are also 
similar to those found by T97, including the velocity overlap 
of the two regions; this overlap is due to the large cloud ve- 
locity dispersion (comparable to the orbital velocity) included 
in the calculations. The emission from the northern and south- 
ern tangent points comprises about half of the total emission, as 
found in T97. The remaining missing flux comes from the east- 
ern and western bridges of smooth emission mostly centered on 
the systemic velocity, as found in POl. The last column of Fig 
4, representing the average over all the 100 realisations, shows 
that the model is stable. That is to say that the main structures 
produced, which match observations, are typical of most reali- 
sations and are not unusual cases. 

Figure 5 shows in more detail the northern tangent region 
for the best matching realisation. The top panel shows good 
overall agreement in the number, distribution and sizes of spots 
as compared with the observations of D99 and T97. Amongst 
these very brightest masers, just as in the observations, there 
are more maser spots on the eastern side of the tangent point 
than the western side. This is explained by the fact that on the 
eastern side there is more background continuum (see Figure 
2a). Note that the middle row of Figure 4 shows an opposite 
behavior on larger scales caused by the presence of the free-free 
absorption bicone. In a given realisation (see inset in FigureS) a 
typical spot spectrum has FWHM of 30 km s ' and FWZI ~ 1 00 
kms ', which is in agreement with the observations presented 
in D99 and T97. The brightest spot shown in Figure 5 is caused 
by the partial overlap of 5 individual clouds. This spot has an 



Parra et al.: A compact starburst ring traced by clumpy OH megamaser emission 



1 1 1 1 1 1 1 1 1 1 1 1 1 
. 1 , . . , 1 , . . , 1 . 


1 1 1 1 1 1 1 1 1 1 1 1 1 
, 1 . . , , 1 , , . . 1 . 


1 1 1 1 1 1 1 1 1 1 1 1 1 
, 1 , . . 1 , 


1 1 1 1 1 1 1 1 1 1 1 1 1 
. 1 , . , , 1 . , , . 1 . 


Ill, . 

■ : 

. 1 , . . , 1 , . . , 1 . 


. , . . . . 1 . . . . , . 
, 1 . . , , 1 , , . . 1 . 


. , . . . . 1 . . . . , . 

' ■ '-h ' ■ " 

. 1 , . . . 1 , , . , 1 . 


., 11 11 1 11 11 1 . 

. \ .... \ .... \ . 



o 
in 



o 
in 



o - 



o 
I 



-50 



50 



-50 



50 -50 
mas 



50 



-50 



50 




-100 100 



100 _-100 
kins 



-100 100 



Fig. 4. The first three columns show sample realisations of the Monte-Carlo simulation. The rightmost column shows the average over 100 
realisations. Intensity and flux are in arbitrary units, normalized to their peaks. Images: Top Row: Greyscale and contours represent the 
velocity integrated OH maser emission convolved with a simulated EVN+MERLIN beam (FWHM » 35 mas). Contours are logarithmically 
spaced from 2"''^ to 2"^'- of the peak. Middle Row: Integrated emission at global VLBI resolution (FWHM = 0.92 mas). Spectra: The solid 
line is the spectrum integrated over the entire source. Dashed and dotted lines represent spectra taken at the southern and northern tangent 
regions respectively. 



LCR of -1200 and a 1667:1665 MHz line ratio p = 28. Both 
values are consistent with the limits LCR > 500 and p > 20 
found by D99 in their brightest spot (see Table 2). 

The position-velocity diagram of the compact maser spots 
shown in the bottom panel of Figure 5 shows a linear gradient 
of ~30±12 kms ' pc ' for the eastern group of spots. This is 
similar to the gradient of 32 km s ' pc ' found by D99 in their 
central group of maser spots (designated N2 in D99). In this 
model realisation no similar gradient was seen in the southern 
tangent region; just as is the case in the observations. This il- 
lustrates the stochastic nature of the mechanism for producing 
velocity gradients as described in §3.6. While our model can re- 
produce the overall gradient seen in the northern tangent region 
the dispersion we find about this gradient is significantly larger 



than that observed by D99. There are several possible reasons 
for this. First it may be that better fits can be found using mod- 
els with larger mean outflow velocities (i.e. larger V; in Table 
1) and smaller random components (i.e. smaller Av and Avj in 
Table 1). Secondly the outflow velocity field we used was very 
crude, simply a change of sign in velocity for clouds on either 
side of the midplane. If clouds are instead gradually acceler- 
ated with increasing z-distance then tighter gradients might be 
produced. An interesting possibility is that such apparent ac- 
celerations could occur if clouds have limited lifetimes before 
their dissipation (see §6.1). In this case clouds emitted from the 
midplane with higher Vz would be observed with systematically 
larger z-distances. Finally our assumption of spherical clouds is 
unlikely to be true in reality. More systematic gradients might 
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Fig. 5. Top: Detail of the northern region of the best matching model 
realisation at global VLBI resolution. The diameter of the circles plot- 
ted are proportional to the maser spot velocity integrated emission. 
For comparison with the results of D99, only spots brighter than 5% 
of the peak are shown. The inset is the spectrum of the brightest feature 
Bottom: Position- Velocity diagram of the brightest maser spots. The 
data suggests the presence of a linear gradient ~30± 12 km s"' pc"' for 
spots on the eastern side of the plot. 

be produced if maser features are instead due to overlapping 
filaments. 

5. Discussion of maser properties 

5. 1. High brightness temperature and maser 
saturation 

Previous studies have argued that compact megamaser features 
imply saturated amplification because of the absence of de- 
tectable seed continuum and the resulting large LCRs (> 500 
in IIIZw35 and > 800 in Arp 220; D99, L98). The high bright- 
ness temperatures of the compact spots, > 2x10'' K in 
IIIZw35 have also been taken as evidence for saturation (D99). 
However, even an LCR of 800 implies maser amplification with 
T of only 6.7, considerably below the typical saturation require- 
ments of astronomical masers, t ~ 12-15 (Elitzur 1992). In our 
model, the brightest features arise from the alignment of « = 5 
maser clouds each with Tc = 1.5 giving a total optical depth 
of 7.5, well below that required for saturation. The assumption 
of unsaturated amplification made in the modelling is therefore 
self consistent even for the strongest features. 

It is also important to remember that saturation depends 
on the angle-averaged intensity J - /x(Q/47r), so even larger 
brightness temperatures and LCRs are possible without satu- 
ration whenever the beaming angles Q. are small. Such small 
beaming angles are natural in the case of a collection of small 
clouds as we have in our model. If there are n spherical masers. 



each of optical depth t^. and diameter a aligned along a length 
L then 

n _ 1 /fl\2 

With a = 0.7 pc, t^. = 1 .5 and L ~13 pc at the ring tangents, 
then the brightest maser spots, for which « = 5, have Q./4n 
of only ~ 2.4x10"^. Hence maser spot brightness temperatures 
many orders of magnitude larger than the present limits would 
still be consistent with unsaturated emission. It is important to 
note that the tight beaming produced by lines of overlapping 
maser clouds is highly anisotropic. The large LCRs and bright- 
ness temperatures achieved in our model would not be observed 
from locations along the ring axis. 

5.2. Maser line ratios 

Compact masers features observed in OH megamasers, display 
large values of p = /leev^Aees- For instance, in Arp 220, L98 
find from their VLBI observations p > 100 in the brightest 
maser spots, in contrast with single dish spectra where p is only 

4.2. In IIIZw35, the 1665 MHz line was detected only toward 
one bright feature at VLBI resolution (D99), yielding p > 20. 
In all other features, this maser emission is below detection. 

The large observed values of p were a major impetus for 
the proposition that the compact emission is different in na- 
ture from the diffuse maser component. However, the observed 
difference in p between these two apparent phases is a natural 
property of clumpy unsaturated masers. If r is the 1667 MHz 
optical depth of a single cloud and R = ti667:ti665, then n over- 
lapping clouds will produce p = exp(nT ■ [1 - 1/^]). From 
an extensive compilation of OH megamaser sources, Henkel & 
Wilson (1990) find that R - 1.9 ± 0.3 for the sample average, 
similar to the ratio of 1.8 if the two transitions have the same 
excitation temperature. For the five-cloud overlap responsible 
for the brightest feature in our simulation, p > 20 for R > 1.7, 
consistent with the observations. The line ratio of the veloc- 
ity integrated emission from our simulated model single dish 
spectra is ~4 which is comparable to the observed value from 
Arecibo spectra of 6.1 (Mirabel & Sanders 1987). 

5.3. Broad maser lines 

In order to fit the observed velocity width of the compact maser 
spots, our model requires clouds with an internal velocity dis- 
persion of 20 km s"' . Although fairly large, this is significantly 
less than if each compact maser were produced in a single 
unsaturated cloud, in which case line narrowing due to the 
high amplification would require cloud velocity widths of ~ 60 
km s ' . In our model, line narrowing due to unsaturated maser 
amplification is largely canceled out by velocity broadening 
from combining many clouds with different velocity centres. 

In bright galactic masers, the brightest spectral features 
are often quite narrow (< 1 kms '). This might be expected 
from a population of individually emitting maser clouds. In 
that case, clouds with small velocity width would give the 
strongest masers because the inverted column density is spread 
over a narrow velocity range so these clouds have the largest 
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peak gains. In contrast, the linewidths of bright megamaser 
spots are typically several tens of kms and up to 150 kms"' 
(Lonsdale 2002). Such large widths cannot arise from a popula- 
tion of narrow-width clouds because they would require unre- 
alistic multi-cloud superpositions. These large widths clearly 
indicate spectral blending of individual clouds with internal 
velocity widths » 1 kms"'. Such blending arises naturally 
from cloud overlaps in a population with a distribution of in- 
ternal cloud properties. Multiple cloud maser emission requires 
overlaps in both space and velocity, therefore narrow-velocity 
clouds, which rarely overlap, are selected against. The over- 
lap of a small number of clouds with large internal velocity 
dispersion, and thus small peak gains, naturally leads to broad 
spectral lines of moderate gain. 

5.4. Application ofttie model to other sources 

Does the model proposed here apply to other sources showing 
both compact and diffuse maser components? The best studied 
OH megamaser source is Arp220, observed at global VLBl res- 
olution (L98) and with EVN and MERLIN separately (Rovilos 
et al. 2003); though there is no published map using the com- 
bined data from both arrays. Arp220 shows maser emission 
from its West and East nuclei, with the latter showing strong 
similarity to IIIZw35: Bright maser emission from two well 
separated regions that contain about a half of the total source 
flux. 

The central velocities of the two bright regions differ by ~ 
100 km s yet their spectral ranges overlap, just as in IIIZw35. 
This suggests that a similar ring model might apply also to the 
Arp220 eastern nucleus; indeed, a rotation of the model im- 
ages presented in Figure 4 shows a striking similarity with the 
Arp220 maps (see Figure 2 in L98). The northern region of 
the eastern nucleus contains more compact structures than the 
southern region, including one very bright spot. This can be 
readily explained in our model, which is stochastic in nature. 
Large velocity gradients over a short distance in the bright re- 
gions are also explained by the overlap of clouds or filaments 
at slightly different centre velocities. In the western nucleus 
of Arp 220 the difference in compactness between the north- 
ern and southern patches of emission are harder to reconcile 
with a rotating ring geometry, although the large LCRs (> 800) 
and large p (> 100) are qualitatively explained by our multiple 
cloud overlap model. 

Another possible candidate source is Mrk273 (Klockner 
et al. 2003). This source, too, shows both compact and 
resolved-out maser emission, and the compact emission occurs 
in two distinct regions with different velocity centres. The ex- 
istence of these two regions can be explained by competition 
between the fall off of background continuum brightness as we 
move away from the projected ring axis versus increased path- 
length through the OH ring. The latter path-length increases 
the maser brightness almost exponentially while the continuum 
strength effects it only linearly. 

Note that in IIIZw35, Mrk273 and Arp220 the brightest 
masers are not coincident with the brightest continuum, but 
this is not strong evidence for saturated emission as is some- 



times claimed. Instead, given the exponential effect of path 
length on brightness in unsaturated masers it is quite feasible 
that the brightest emission occurs in regions where the contin- 
uum emission is very weak or presently undetectable. 

6. Physical considerations 

6. 1. Cloud mass and confinement 

Here we discuss the physical properties of our OH maser 
clouds. The maximum density of a OH maser emitting gas is of 
order «(H2) = 10^ cm"^. Higher densities will thermalise the 
energy levels and quench the maser (Elitzur 1992). Assuming 
this maximum density, a OH abundance of 10"^ and excitation 
temperature lOK, equation 9.12 in Elitzur (1992) gives a min- 
imum cloud diameter of 0.02pc. Assuming instead that clouds 
have diameters equal to the observational upper limit of 0.7pc 
then the same equation gives a minimum mean hydrogen den- 
sity of 3.1xlO-'cm"''. This gives an upper limit for cloud mass 
of 24M0. Note however that this upper limit is critically depen- 
dant on the assumed OH abundance. 

Are the OH maser clouds confined? It seems that gravi- 
tational confinement can be rejected. Given their internal ve- 
locity dispersion of 20km s"', to be gravitationally bound the 
clouds would require a virial mass of ~2xlO'*Mo. This is 
much larger than the upper limit estimated from OH proper- 
ties. Such a cloud mass would also give a total mass in clouds 
of ~2xlO^M0, which is larger than the dynamical mass internal 
to the maser ring of VxIO^Mq derived from the maser kinemat- 
ics (POl). 

Could the clouds be pressure confined by the ionised Inter- 
Cloud Medium (ICM)? Using the inferred optical depth of 
the free-free absorbing gas and assuming it is uniformly dis- 
tributed, the electron number density in the ICM is ~10'' cm"^. 
If the temperature of the ICM is 10* K, then the pressure be- 
comes PicM~10^ K cm"^. In contrast, from the internal tur- 
bulent velocity of the clouds their inferred dynamical pres- 
sure is of order 3x10** K cm therefore they cannot be pres- 
sure confined. Finally magnetic confinement would require a 
magnetic field of ~10mG. This is feasible given the range of 
magnetic field strengths observed towards galactic OH masers 
(Reid & Moran (1981), Fiebig & Guesten (1989)) if the OH 
masers occurred in gas of the highest possible density (i.e near 
n(H2) = 10^ cm"^). For comparison in IIIZw35 Zeeman split- 
ting observations by Killeen et al. (1996) estimate an upper 
limit on the non-uniform magnetic field to be 5mG, however 
this limit is very model dependant. 

Although magnetic confinement is possible given the ob- 
servational limits, it is interesting to consider models in which 
the clouds are freely expanding. Assuming the largest possible 
cloud size (0.7pc) and dividing by the cloud internal velocity 
dispersion gives a characteristic cloud lifetime of ~3.4xl0'* yr. 
Given their outflow velocity of 60km s ' in the z direction, such 
clouds would reach a height of ~2. 1 pc which is roughly consis- 
tent with the observed A///2=3pc (see Table 1). An additional 
advantage of this model is that that the ionised free-free ab- 
sorbing gas required by our model might naturally be generated 
from the dissipated clouds. If the clouds have uniform density 
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and consequently a mass of ~24Mq then the total kinetic power 
being injected into all clouds would be about TxlO-'^ erg s"' 
which is a very small fraction of the mechanical power avail- 
able from supernova explosions (i.e. 2x10'*-' erg s"') assuming 
a supernova rate of 0.8yr"' (POl). Using the same cloud mass, 
the mass loss rate in the outflowing clouds is O.SMoyr"' which 
is much smaller than the estimated SFR (19M0yr"', see POl) 

6.2. Constraints on central AGN and Black Hole 

Our axisymmetric model (Figure 2) suggests a central point to 
both the ring and the ionised outflow and it is natural to ask 
whether an AGN or black hole exists at this point. We should 
first remark that within the OH ring radius, the geometry of our 
model is observationally ill-defined and so the cone of free- 
free absorbing gas shown in figure 2 may not in fact extend all 
the way inwards to a central point. If as we argue in §6.1, the 
free-free absorbing gas is the remains of dissipated molecular 
clouds, then the free-free cone will in fact be truncated at a 
height comparable to the OH ring height. 

Observationally there is conflicting evidence for the pres- 
ence of an AGN in lllZw35. The source lies on the well known 
FlR-radio correlation for starbursts (see POl) and so has no ra- 
dio excess which might accompany a strong AGN. Also high 
resolution radio observations (see POl) do not show any com- 
pact radio features at our inferred centre. However such a ra- 
dio core might be free-free absorbed at 1.6GHz so it would be 
interesting to conduct higher frequency VLBl observations to 
check this possibility. We have not been able to find any X-ray 
observations which imply an AGN. 

In contrast to the above C90 argued that the near-lR 
colours of lllZw35 were consistent with it being a mixed star- 
burst/ AGN. Furthermore from their optical spectroscopy C90 
classified lllZw35 as a borderline L1NER/Seyfert2 nucleus 
based on classical (semi-empirical) line ratio diagnostics. A 
similar analysis by Baan et al. (1998) using different opti- 
cal spectroscopic observations but a similar diagnostic scheme 
classified lllZw35 as a LINER. Whether galaxies classified as 
LlNERs are primarily AGN or starburst powered is presently 
unclear Using the physically based line diagnostic scheme of 
Kewley et al. (2001) many LlNERs are reclassified as high 
metallicity starbursts (Corbett et al. 2003). However, taking 
the data for lllZw35 from Baan et al. (1998) and applying the 
Kewley et al. (2001) diagnostic scheme we find that lllZw35 
still falls in the AGN region. The data is however consistent 
with a mixed AGN/starburst with the latter still providing most 
of the luminosity in the diagnostic optical lines. As noted by 
Corbett et al. (2003). estimating the relative bolometric AGN 
contribution in starbursts purely from optical emission lines is 
very difficult because of the large optical obscuration in these 
objects. 

Is there any dynamical evidence for a central black hole? 
From the observed maser dynamics the derived enclosed mass 
within the OH maser ring radius of 22pc is TxIO^Mq . The re- 
sulting mass density of I56M0 pc"^ can easily be achieved by 
stars or gas within the nuclear region and so there is no neces- 
sity for a central point mass. However a moderate mass black 



hole is not ruled out either The total enclosed mass in IIIZw35 
could for instance be contributed equally from a black hole like 
the one in our galaxy (3.5xlO''M0 , Eckart & Genzel (1997)) 
and a distributed mass of stars. We note that even in the ex- 
treme case that such a black hole is being maximally fed, then 
the maximum luminosity if radiating at 10% of the Eddington 
limit would be 2xlO'"L0 which is over ten times smaller than 
the observed FIR luminosity. We conclude that although optical 
spectroscopy shows that an AGN might be present in IIIZw35 
it must be energetically insignificant compared to the starburst 
activity. 

6.3. Nature of the OH and continuum rings 

We argue in the previous section that the bulk of the bolomet- 
ric emission in IIIZw35 is powered by starburst activity rather 
than an AGN. Certainly for the radio emission its consistency 
with the radio-FIR correlation, its observed brightness temper- 
ature (POl) and its ring-like morphology argues for an origin in 
a starburst ring. Such a starburst ring can also provide the IR 
photons to pump OH and therefore also explain the ring mor- 
phology of the maser emission. 

A remarkable property of the OH ring is its relative narrow- 
ness (see §3.3 and Table 1), which is only 12% of the radius. 
The fact that OH absorption is seen outside the ring and perhaps 
in its central hole (POl) suggests that this geometry is not de- 
fined merely by OH abundance. The most direct interpretation 
of the narrow maser ring is that it is defined by the range of radii 
over which star formation is presently occurring. Such ring can 
be identified as a scaled down version of the sub-kiloparsec 
nuclear starburst rings known to exist in many starburst galax- 
ies (see Knapen et al. (2004) and references therein). In addi- 
tion, numerical simulations of starburst-generated rings or tori 
have been made (Wada & Norman (2003), Wada & Tomisaka 
(2005)) supporting the existence of circumnuclear filamentary 
structures. The presence of radio continuum emission at larger 
radii in IIIZw35 (see Figure 1) may indicate the remnant of 
past star formation, in turn suggesting an inwardly propagating 
ring of star-formation. The fact that the OH masing clouds in 
the ring are not gravitationally bound to the nucleus (see §6.1) 
suggests that emission at any given radius must be short lived. 
Supernova explosions in a ring may compress dense gas on the 
inside of the ring causing the ring to propagate inwards. 

Alternatively to the above picture, star formation in 
IIIZw35 may exist over a wide range of radii, rather than just 
in a ring, but the required conditions for population inversion 
might only apply over a narrow range of radii. The inversion of 
OH main lines is critically dependant on the local IR spectrum 
(EUtzur (1978), Cohen et al. (1988), Moore et al. (1988), Kegel 
et al. (1999)), which in turn depends on the local SFR rate and 
Initial Mass Function. 

It is interesting to compare the size and structure of our OH 
maser ring with other OH maser observations of circumnuclear 
features. In Mrk231 (Klockner et al. 2003) an OH torus with 
radius from 30pc to lOOpc has been inferred. It follows that 
the inner edge of that torus is comparable in size to the ring 
in IIIZw35 but it extends over a larger radius and appears to 
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be much thicker. We argue in §5.4 that the structures seen at 
the eastern nucleus of Arp220 (Rovilos et al. 2003) could also 
be explained by a thin ring with a radius of ~30pc being very 
similar to IIIZw35. It appears from the observations that a wide 
variety of circumnulear OH maser structures can exist. 

The OH torus found in Mrk 231 (Klockner et al. 2003) has 
been identified with the obscuring torus required in AGN uni- 
fied schemes, although the inner radii of such structures are 
expected to be much smaller (close to the dust sublimation ra- 
dius which is sub-parsec in weak Seyferts or strong LINERS). 
In contrast it does not seem that the OH ring we see in lllZw35 
can directly contribute much to any such obscuration. From 
the direction of observation, the ring does not cover the central 
point where any AGN might be expected to reside. Furthermore 
the naiTow ring covers only a small solid angle and even if it 
was observed edge-on, the obscuration would still be small be- 
cause given the derived cloud column densities (§6.1) the ex- 
pected obscuration is only A,,~l assuming a standard gas to 
dust ratio. However, it cannot be ruled out that the structure 
seen in OH is the outer part of a starburst supported disk which 
extends down to small radii where it might have a much larger 
geometrical and optical thickness. 

7. Conclusions and future work 

Observations of distinct regions of compact and apparently dif- 
fuse maser emission in OH megamasers have, in the past, been 
used to argue for two physical phases of OH masers (see L98, 
D99). In contrast, the mechanism discussed in this paper ex- 
plains both types of maser structures using a single phase of 
low opacity clouds within a thin circumnuclear ring. The prop- 
erties of the clouds are similar to those assumed in the standard 
model (Baan 1989) but our model explicitly considers the sta- 
tistical effects of rare multi-cloud overlaps. Compared to oth- 
ers, our model is therefore conservative in that the physical 
properties of the standard model are preserved. 

Despite the simplicity of our model we find that most of the 
observational features of IIIZw35 can be reproduced. The fact 
that the maser amplifying medium is composed of clouds is 
found essential to explain the range of maser brightness around 
the ring at low resolution. The same clouds also explain the 
bright maser spots seen only at the ring tangents in terms of 
multiple cloud overlaps in both space and velocity. The model 
is able to reproduce the LCRs in the bright spots and diffuse re- 
gions, and the large value of the 1667MHz: 1665MHz line ratio 
in the compact spots. Finally it can explain the spectra of both 
the compact spots and the apparently diffuse areas of emission. 
The fact that the OH clouds are outflowing from the ring mid- 
plane explains the large velocity gradients observed amongst 
the compact maser spots. We find that the OH masing ring is 
relatively narrow in radius which could be explained either in 
terms of a narrow circumnuclear ring of star formation, or due 
to the strong sensitivity of maser pumping to physical condi- 
tions which vary gradually with radius. The ring we find is nar- 
rower but qualitatively similar to those that have been produced 
in numerical simulations (Wada & Norman (2002), Wada & 
Tomisaka (2005)). 



Our model can be improved in many ways. It assumes 
all clouds are identical, when in reality a spectrum of cloud 
sizes and opacities is expected. Additionally it presently as- 
sumes spherical clouds whereas computer simulations (Wada 
& Norman (2002)) show that filamentary structures are com- 
mon in starburst nuclei. These latter structures may also show 
outward velocity gradients due to acceleration by radiation 
pressure; including such filaments will give rise to higher ve- 
locity correlation between masers spots, as seen in the observa- 
tions. The next step in simulations of megamaser sources might 
be to take the output of numerical dynamical simulations and 
calculate the expected OH megamaser emission. 

One line of future work involves taking the estimates of 
cloud opacity, size, internal velocity dispersion and 1667:1665 
MHz line ratio found in IIIZw35 and deriving the implications 
for gas physical conditions and pumping. Another important 
area to investigate is the general one of properties of maser 
emission in a cloudy medium, especially the generalisation to 
media with a range of cloud opacities. While much work has 
been put into understanding masers in simple geometries such 
as spheres or filaments relatively little work has been done for 
random media. One exception is Sobolev et al. (2003) who car- 
ried out numerical simulations of masers in continuous media 
with a Kolmogorov spectrum of local opacity. Multi -phase me- 
dia are common in astrophysics and the case we have studied 
here of masers arising in discrete clouds is another limit which 
requires careful attention. Observationally separate regions of 
compact and apparently diffuse maser emission are often found 
in other species of galactic and circumstellar masers. Examples 
include galactic methanol (Minier et al. 2000) and circumstel- 
lar SiO masers (Yi et al. 2005). These core-halo structures are 
often interpreted in terms of single clouds which have saturated 
outer layers and unsaturated cores. However, such structures 
might also be explained by invoking the presence of a clumpy 
medium. 

Appendix A: Maser gain in ciumpy media 

Consider a population of identical unsaturated maser clouds 
amplifying background continuum. At a given velocity let 
be the mean number of clouds which are encountered along a 
given LOS. If the clouds are of non-zero size and cannot inter- 
penetrate there is also a maximum number of clouds A^max, set 
by the number than can be fitted along a LOS. If the maser op- 
tical depth of each cloud is Tc then the mean value of the gain 
is 

A'mux 

A-=0 

where (A:) is the probability of having k overlaps. If the 

cloud volume filling factor is small such that A^max >> N then 
Pfi f^^ik) can be replaced by the Poisson distribution /'^(^) - 
fjkg-N and the upper limit on the sum tends to infinity. In 
this case the mean of the gain is given by 

G„ean= ^ P^(^) ■ e*^^' = cxp [^(e^' - 1)] (A.2) 

k=Q 
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as found by Natta & Panagia (1984). However we find that 
when Tc > 1 the above expression enormously overestimates 
the average gains in our Monte-Carlo simulations; even for 
a negligible cloud filling factor The reason is that in the 
above series although the Poisson probabilities Pj^{k) decreases 
rapidly for k > N, for Tc > 1 this decrease is more than offset by 
the increase of the amplification factors e*'^'^ . The result is that 
the series summation becomes dominated by very rare events 
with k » N. In contrast for a typical Monte-Carlo simulation 
such an extreme multi-cloud overlap event is very unlikely to 
occur even once within any of the independent lines of sight 
within an interferometer beam. 

An analytic expression for Gmean which is in closer agree- 
ment with our Monte-Carlo simulations can be constructed in 
the following way: If the interferometer beam is M times larger 
than a cloud area, each beam can be assumed to contain M in- 
dependent lines of sight. Consequently, the series in equation 
A. 2 can be truncated at a finite A^max to exclude improbable 
events that are unlikely to occur even once within the beam. 
More accurately, A^max can be defined such that the maximum 
number of overlaps within a region of M lines of sight is less 
than Atnax in 50% of realisations. This gives: 

Gmean W ■ e*"^ A^max = (2^^) (A-3) 

where P^} is the inverse Poisson cumulative distribution. Using 

N ° 

Monte-Carlo simulations (see Figure A.l) this equation is 
found to give a good approximation to the mean gain over a 
wide range of N. Using the above definition of Amax the me- 
dian value of the peak gain within a region of M independent 
lines of sight is given by 

Gmax - eXp(Aniax ' Tc) (A.4) 

An important point, which is evident from Figure A.l, is 
that the effective opacity Tea = In(Gniean) shown on the scale 
to the right, is much larger than the mean opacity in that region 
Tmean = Nt^ (showu at the scale on the top). The two quantities 
differ because of the high k tail of the f ^ distribution. Figure 
A. 1 also demonstrates another remarkable property of masers 
in a clumpy medium, that the effective opacity Tefj, is not lin- 
early dependent on N but falls well below that prediction as 
N and Tmean increase. This effect arises because as N increases 
the variance about the mean in the number of clouds per line 
of sight decreases which in turn reduces the contribution from 
rare many-cloud overlaps. As described in §3.4 this nonlinear 
relationship between effective opacity and .^V for a clumpy me- 
dia allows us to understand the otherwise puzzling variation in 
LCR around the maser ring. 
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Fig. A.l. Variation of the mean gain G with the mean number N of 
overlapping maser clouds when the individual cloud optical depth is 
Tc = 1.5. The crosses show the results of a numerical Monte-Carlo 
calculation that averages over M = 50 lines of sight. The solid line is 
the result of analytic calculation using equation A. 3. The dashed line 
plots the result of an infinite-series summation (Equation A. 2). The 
top axis shows the distribution mean optical depth Tmem = Nt^, the 
right axis its effective optical depth T^ff = lnG„,ean- 
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